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ABSTRACT 

This paper describes a new grid for the mass-radius relation of 3-layer exoplan- 
ets within the mass range of 0.1 through 100 M^. The 3 layers are: Fe (e-phase of 
iron), MgSiOs (including both the perovskite phase, post-perovskite phase, and 
its dissociation at ultra-high pressures), and H2O (including Ices Ih, III, V, VI, 
VII, X, and the superionic phase along the melting curve). We discuss the cur- 
rent state of knowledge about the equations of state (EOS) that influence these 
calculations and the improvements used in the new grid. For the 2-layer model, 
we demonstrate the utility of contours on the mass-radius diagrams. Given the 
mass and radius input, these contours can be used to quickly determine the im- 
portant physical properties of a planet including its pO (central pressure), pl/pO 
(core- mantle boundary pressure over central pressure), CMF (core mass fraction) 
or CRF (core radius fraction). For the 3-layer model, a curve segment on the 
ternary diagram represents all possible relative mass proportions of the 3 layers 
for a given mass-radius input. These ternary diagrams are tabulated into Table [3] 
with the intent to make comparison to observations easier. How the presence of 
Fe in the mantle affects the mass-radius relations is also discussed in a separate 
section. A dynamic and interactive tool to characterize and illustrate the inte- 
rior structure of exoplanets built upon models in this paper is available on the 
website: ,http : / / www . cf a . harvard . edu/ ~ Izeng, 

Subject headings: exoplanet, Super-Earth, interior structure modeling, Kepler 
Mission 



1. Introduction 



The transit method of exoplanet discovery has produced a small, but well constrained, 
sample of exoplanets that are unambiguously solid in terms of interior bulk composition. We 
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call solid planets the ones that possess no H and He envelopes and/or atmospheres, i.e. their 
bulk radius is determined by elements (and their mineral phases) heav ier than H and He . 



Such solid exoplan ets are Keple r-lOb (iBatalha et al.ll201ll ). CoRoT-7b f lQueloz et al.l 12009 
Hatzes et al. 2011 ). Kepler-36b ( Carter et al. 20121 ) as well as - most likely: Kepler-20b,e,f; 



Kepler-18b, and 55 Cnc e in which the solid material could include high-pressure water ice 
(see references to Fig. I3]in l4.2p . 

There is an increased interest to compare their observed parameters to current models 
of interior planetar y structure. The models, and their use of approximations and EOS, have 
evolv ed since 2005 ( Valencia et al.|[2006 : Fortney et al.lboOT : Grasset et al. 2009 : Seager et al. 
20071 ). mainly because the more massive solid exoplanets (called Super-Earths) have interior 
pressures that are far in excess of Earth's model, bringing about co rresponding gaps in our 
knowledge of mineral phases and their EOS (see a recent review by lSotin et al.ll2010[ ). Here 
we compute a new grid of models in order to aid current comparisons to observed exoplanets 
on the mass-radius diagram. As in previous such grids, we assume the main constituents 
inside the planets to be differentiated and model them as layers in one dimension. 

The first part of this paper aims to solve the 2-layer exoplanet model. The 2-layer model 
reveals the underlying physics of planetary interior more intuitively, for which we consider 3 
scenarios: Fe-MgSiOs or Fe-H20 or MgSi03-H20 planet. 

The current observations generally measure the radius of an exoplanet through transits 
and the mass through Doppler shift measurement of the host star. For each assumption of 
core and mantle compositions, given the mass and radius input, the 2-layer exoplanet model 
can be solved uniquely. It is a unique solution of radial dependence of interior pressure and 
density. As a result, all the characteristic physical quantities, such as the pressure at the 
center (pO), the pressure at core-mantle boundary (pi), the core mass fraction (CMF), and 
the core radius fraction (CRF) naturally fall out from this model. These quantities can be 
quickly determined by invoking the mass-radius contours. 

The next part of this paper (Fig. [3]) compares some known exoplanets to the mass- 
radius curves of 6 different 2-layer exoplanet models: pure-Fe, 50%-Fe & 50%-MgSiO3, pure 
MgSiOa, 50%-MgSiO3 & 50%-H2O, 25%-MgSi03 & 75%-H20, and pure H2O. These per- 
centages are in mass fractions. The data of these six curves are available in Table [1] 

Up to now, a standard assumption has been that the planet interior is fully differentiated 
into layers: all the Fe is in the core and all the MgSiOa is in the mantle. In section 14. 3[ 
we will change this assumption and discuss how the presence of Fe in the mantle affects the 
mass-radius relation. 



The final part of this paper calculates the 3-layer differentiated exoplanet model. Given 
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the mass and radius input, the solution for the 3-layer model is non-unique (degenerate), 
thus a curve on the ternary diagram is needed to represent the set of all solutions (see Fig. |5]). 
This curve can be obtained by solving differential equations with iterative methods. The 
ensemble of solutions is tabulated (Table |3]), from which users may interpolate to determine 
planet composition in 3-layer model. A dynamic and interactive tool to characterize and 
illustrate the interior structure of exoplanets built upon Table [3] and other models in this 
paper is available on the website http: //www, cf a .harvard . edu/~lzen g, 

The methods described in this paper can be used to fast characterize the interior struc- 
ture of exoplanets. 



2. Method 



Spherical symmetry is assumed in all the models. The interior of a planet is assumed to 
be fully differentiated into layers in the first part of the paper. The 2-layer model consists of 
a core and a mantle. The 3-layer model consists of a core, a mantle and another layer on top 
of the mantle. The interior structure is determined by solving the following two differential 
equations: 



dr 


1 


dm 




dp 


Gm 


dm 


47rr^ 



(1) 



(2) 



The two equations are similar to the ones in IZeng fc Seagen (120081 ). However, contrary 
to the common choice of radius r as the independent variable, the interior mass m is chosen, 
which is the total mass inside shell radius r, as the independent variable. So the solution is 
given as r(m) (interior radius r as a dependent function of interior mass m), p(m) (pressure 
as a dependent function of interior mass m), and p(m) (density as a dependent function of 
interior mass m). 



The two differential equations are solved with the EOS of the constituent materials as 
inputs: 



P = P{P,T) (3) 
The EOS is a material property, which describes the material's density as a function 



-4- 



of pressure and temperature. The EOS could be obtained both theoretically and exper- 
imentally. Theoretically, the EOS could be calculated by Quantum-Mechanical Molecu- 
lar D ynamics Ab Initio Sirnulation such as the VASP (Vienna Ab initio Simulati on Pack- 
age) JKresse fc Hafnerlll993l . Il994l : iKresse fc Furthmullerlll996l : IPrench et aPboOQi ). Exper- 
imentally, the EOS could be determined by high-pressure compression experiment such as 
the DAC (Diamond Anvil Cell) experi ment, or shock wave experiment like the implosion 
experiment by the Sandia Z-machine (lYu fc JacobsenI 120111). The temperature effect on 
density is secondary compared to the pressure effect (I Valencia et al.l l2006l . l2007bl ) . There- 
fore, we can safely ignore the temperature dependence of those higher density materials 
(Fe and MgSiOs) for which the temperature effect is weaker, or we can implicitly include a 
pre-assumed pressure-temperature (p-T) relation (for H2O it is the melting curve) so as to 
reduce the EOS to a simpler single- variable form: 



P = pip) 



(4) 



To solve the set of equations mentioned above, appropriate boundaries conditions are 



given as: 



• pO: the pressure at the center of the planet 

• pi: the pressure at the first layer interface (the core-mantle boundary) 

• p2: the pressure at the second layer interface (only needed for 3-layer model) 

• Vsurface- the prcssure at the surface of the planet (set to 1 bar (lO^Pa)) 

3. EOS of Fe, MgSiOa and H2O 

The 3 layers that we consider for the planet interior are Fe, MgSiOs and H2O. Their 
detailed EOS are described as follows: 

3.1. Fe 

We model the core of a solid exoplanet after the Earth's iron core, except that in our 
model we ignore the presence of all other elements such as Nickel ( Ni) and light e l ement s 
such as Sulfur (S) and Oxygen (O) in the core. As pointed out by [Valencia et al.l ( I2OIOI ). 
above 100 GPa, the iron is mostly in the hexagonal closed packed e phase. Therefore, we 
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use the Fe- EOS bv ISeager et all (120071 ). Below 2.09 * 10^ GPa, it is a Vi net fIVinet et al. 
987l. Il989l ) formula fit to the experimental data of e-iron at p< 330 GPa by lAnderson et al.l 



(12001 



Above 2.09*10^ GPa, it makes smooth transition to the Thomas- Fermi- Dirac (TFD) 



EOS (jSeager et al.l 120071 ). A smooth transition is assumed because there is no experimental 



data available in this ultrahigh-pressure regime. 

The central pressure could reach 250 TPa (terapascal, i.e., 10^^ Pa) in the most massive 
planet considered in this paper. However, the EOS of Fe above 400 GPa is beyond the 
current reach of experiment and thus largely unknown. Therefore, our best approximation 
here is to extend the currently available e-iron EOS to higher pressures and connect it to the 
TFD EOS. 



The EOS of Fe is shown in Fig. [T]as the upper curve (red curve). 



3.2. MgSiOa 



We model the silicate layer of a solid exoplanet using the Earth's mantle as a proxy. The 
FeO-free Earth's mantle with Mg/Si=1.07 would consist of mainly enstatite (MgSiOs) or its 
high-pressure polymorphs and, depending upon pressure, small amou nts of eithe r forsterite 
and its high-pressure polymorphs (Mg2Si04) or periclase (MgO) (e.g. 



Binall2003h 



The olivine polymorphs as well as lower-pressure enstatite and majorite (MgSiOa with 

garnet structure), are not stable above 27 GPa. At higher pressures, the syste m would consist 

o f MgSiO.s-perov skite (pv) and periclase or their higher-pressure polymorphs ( IStixrude fc Lithgow-Bertellori 



2011 



Bina 



20031 ). Given the high pressures at the H20-silicate boundary usually exceeding 



27 GPa, we can safely ignore olivine and lower-pressure pyroxene polymorphs. For the sake 
of simplicity, we also ignore periclase, which would contribute only 7 at.% to the silicate 
mantle mineralogical composition. There are also small amounts of other elements such 



as Al uminum (Al), Calcium (Ca), and Sodium (Na) present in Earth's mantle (ISotin et al. 



20071 ). For simplicity, we neglect them and thus the phases containing them are not included 
in our model. We also do not consider SiC because carbon-rich planets might form under 
very rare circumstances, and are probably not common. 

Some fraction of Fe can be incorporated into the minerals of the silicate mantle which 
could then have the general formula as (Mg,Fe)Si03. For now, we simply assume all the Fe 
is in the core and all the Mg is in the mantle in the form of MgSiOs-perovskite and/or its 
high-pressure polymorphs. So the planet is fully differentiated. In a later section, we will 
discuss how the addition of Fe to the mantle can affect mass-radius relation and compare 
the differences between differentiated and undifferentiated as well as reduced and oxidized 
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p (pressure, Pa) 



Fig. 1. — P — p EOS of Fe (e-phase of iron, red curve), MgSiOa (perovskite phase, post- 
perovskite phase and its high-pressure derivatives, orange curve), and H2O (Ice Ih, Ice III, 
Ice V, Ice VI, Ice VII, Ice X, and superionic phase along its melting curve (solid-liquid phase 
boundary)) 
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planets in Section 14.31 



We first consider the perovskite (pv) and post-perovskite (ppv) phases of pure MgSiOs. 
MgSiOs perovskite (pv) is believed to be the major constituent of the Earth mantle. It 
makes transition into the post-perovskite (ppv) phase a t roughly 120 GPa and 2500 K (cor- 
responding to a depth of 2600 kilo meters in Earth! (jHirosd 120101 ). The ppv phase was 



pv phase (ICaracas fc Cohen 



2008 



Hirose 



discovered experimentally in 2004 fbvlMurakami et al.ll2004f) and was also theoretically pre- 
dicted in t he same year (by Qganov fc Oiidl2004r) . The ppv is about 1.5% denser than the 



20101). This 1.5% density jump resulting from 



the pv-to-ppv phase transition can be clearly seen as the first density jump of the MgSiOs 
EOS curve shown in F i g. [H Both the MgSiOs pv EOS and MgSiOa ppv EOS are taken 
from ICaracas fc CohenI (120081 ). The t r ansiti on pressure is determined to be 122 GPa for 



pure MgSiOs according to ISpera et al.l (120061 ) 



Beyond 0.90 TPa, MgSiOs ppv u ndergoes a two-stage dissociation process predicted 
from the first-principle calculations by lUmemoto fc WentzcovitchI (120111 ). MgSiOs ppv first 
dissociates into CsCl-type MgO and P2ic-type MgSi205 at the pressure of 0.90 TPa and later 

into CsCl-type MgO and Fe2P-type Si02 at pressures higher than 2.10 TPa. The EOS of 

CsCl- ty pe MgO, P2ic-ty pe MgSi205, and Fe2P-type Si02 are adopted from lUmemoto fc Wentzcovitch 
(I2OIII ): IWu et al.l (l201l[ ). Therefore, there are two density jumps at the dissociation pres- 
sures of 0.90 TPa and 2.10 TPa. The first one can be seen clearly in Fig. [H The second one 
cannot be seen in Fig. [T] since it is too small, but it surely exists. 



Since lUmemoto fc WentzcovitchI (1201 II ) ' s EOS calcul a tion is until 4.90 TPa, beyond 
4.90 TPa, a modified version of the EOS by ISeager et al.l ( 120071 ) is used to smoothly con- 
nect to the TFD EOS. TFD EOS assumes electrons in a slowly varying potential with a 
density- dependent correlation energy term that describes the interactions among electrons. 
It is therefore insensitive to any crystal structure or arrangements of atoms and it becomes 
asymptotically more accurate at higher pressure. Thus, the TFD EOS of MgSiOs would 
be no different from the TFD EOS of MgO plus Si02 as long as the types and numbers of 
atoms in the calculation are the same. So it is safe to use the TFD EOS of MgSiOs as an 
approximation of the EOS of MgO and Si02 mixture beyond 4.90 TPa here. 



Seager et al.l ( 120071 ) 's EOS is calculated from the method of ISalpeter fc ZapolskyI (Il967| ) 
above 1.35 * 10^ GPa. Below 1.35 * 10^ GPa, it is a smooth connection to TFD EOS 



from the fourth-order Birch-Murbaghan Equation of State (BME) (see iBirchI Il947l : iPoirier 
20001) fit to the param eters of MgSiOs pv obtained by Ab initio lattice dynamics simulation 
of iKarki et aP tood ). 



Seager et al, 



's EOS is slightly modified to avoid any artificial density jump when con- 
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nected with lUmemoto fc Wentzcovitchl's EO S at 4.90 TPa. At 4.90 TPa, the ratio of the 
density p between lUmemoto &: Wentzcovitchl's EOS and lSeager et aLr s EOS is 1.04437. This 
ratio is multiphed to the original ISeager et aLr s EOS density p to produce the actual EOS 
used in our calculation for p> 4.90 TPa. A smooth transition is assumed because no exper- 
imental data is available in this ultrahigh-pressure regime. This assumption does not affect 
our low or medium mass planet models, since only the most massive planets in our model 
could reach this ultrahigh pressure in their MgSiOa part. 

The EOS of MgSiOs is shown in Fig. [T]as the middle curve (orange curve). 



3.3. HoO 



The top layer of a planet could consist of various phases of H2O. Since H2O has a complex 
phase diagram, and it also has a stronger temperature dependence, thus the temperature 
effect cannot be ignored. Instead, we foll ow the solid ph ases along the melting curve (solid- 
liquid phase boundary on the p-T plot by I Chaplin! ( 120121 )). Along the melting curve, the H2O 
undergoes several phase transitions. Initially, it is Ice Ih at low pressure, the n subsequently 



transforms into Ice III, Ice V, Ice VI, Ice VII, Ice X, and superioni c phase (IChaplinI 12012 
Choukroun fc Grassetl [2OO7I : IPunaeva et aPboiol : IPrench et al]l2009[ l. 



3.3.1. Chaplin's EOS 



The solid form of water has very complex phases in the low-pressure and low-temperature 
regime. These phases are well determined by experimen ts. Here we adopt the Chaplin's EOS 
for Ice Ih, Ice III, Ice V, and Ice VI below 2.216 GPa (see lChaplinll2012l : IChoukroun fc Grasset 
20071 : iDunaeva et al.ll2010l ). Along the melting curve (the solid-liquid boundary on the p-T 
diagram), the soli d form of water first exists as Ice Ih (Hexagonal Ice) from ambient pressure 
up to 209.5 MPa Jchoukroun fc Grassetl l2007h. At the triple point of 209.5 MPa and 251.15 



K (I Choukroun fc Grasset 



2007 



Robinson et al.lll996l ). it transforms into Ice HI (Ice-three), 



whose unit cell forms tetragonal crystals. Ice HI exists up to 355.0 MPa and transforms 
int o a higher-pressure form Ice V (Ice-five) at the triple point of 355.0 MPa and 256.43 
K (jchoukroun fc Grassetl 120071 ). I ce V's unit cell forms monoc linic crystals. At the triple 
point of 618.4 MPa and 272.73 K (jChoukroun fc Grassetl 120071 ) . Ice V transforms into yet 
another higher-pressure form Ice VI (Ice-six). Ice VI's unit cell forms tetragonal crystals. 
A single molecule in Ice VI crystal is bonded to four other water molecules. Then at the 
triple point of 2.216 GPa and 355 K (jPaucik &: Dooleyll201ll ). Ice VI transforms into Ice VII 
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(Ice-seven). Ice VII has a cubic crystal structure. Ic e VII eventually transf orms into Ice X 
(Ice-ten) at the triple point of 47 GPa and 1000 K (iGoncharov et al.ll2005l ). In Ice X, the 
protons are equally spaced and bond ed between the oxygen ato ms, where the oxygen atoms 
are in a body-centered cubic lattice ( iHirsch &: Holzapfellll984j ). The EO S of Ice X and Ice 
VII a re very similar. For Ice VII (above 2.216 GPa), we switch to the (IFrank. Fei. fc Hu 
2004h 's EOS (FFH2004). 



3.3.2. FFH2004's EOS 

We adopt FFH2004's EOS of Ice VII for 2.216 GPa<p< 37.4 GPa. This EOS is obtained 
using the Mao-Bell type diamond anvil cell with an external Mo-wire resistance heater. Gold 
and water are put into the sample chamber and compressed. The diffraction pattern of both 
H2O and gold are measured by the Energy-Dispersive X-ray Diffraction (EDXD) technique 
at the Brookhaven National Synchrotron Light Source. The gold here is used as an internal 
pressure calibrant. The disappearance of the diffraction pattern of the crystal Ice VII is 
used as the indicator for the solid-liquid boundary (melting curve). The melting curve for 
Ice VII is determined accurately from 3 GPa to 60 GPa and fit to a Simon equation. The 
molar density (p in mol/cm?) of Ice VII as a function of pressure (p in GPa) is given by the 
following formula in FFH2004: 



p{mol/cm^) = 0.0805 + 0.0229 * (1 - exp 



'0.0743*p 



) + 0.1573 *(1 - exp 



-0.0061*p 



) (5) 



We use Eq. [5]to calculate p from 2.216 GPa up to 37.4 G Pa. The upper limit 37.4 GPa is 

deter mined by the intersection between FFH2004's EOS and (IFrench. Mattsson. Nettelmann. fc Redmer 
2o3)'s EOS (FMNR2009). 



3.3.3. FMNR2009's EOS 



FMNR2009's EOS is used for Ice VII, Ice X and superionic phase of H2O for 37.4 
GPa<p< 8.893 TPa. This EOS is determined by Quantum Molecular Dynamics Simulations 
using the Vienna Ah Initio Simulation Package (VASP). The simulation is based on finite 
temperature density-functional t heory (DFT) for the electronic structure and treating the 
ions as classical particles. Most of French et al.r s simulations consider 54 H2O molecules in a 
canonical ensemble, with the standard VASP PAW potentials, the 900 eV plane-wave cutoff, 
and the evaluation of the electronic states at the F point considered, for the 3 independent 
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variables: temperature (T), volume (V), and particle number (N). The simulation results 
are the thermal EOS p(T,V,N), and the caloric EOS U(T,V,N). The data are tabulated in 
FMNR2009. 

In order to approximate density p along the melting curve, from Table V in FMNR2009, 
we take 1 data point from Ice VII phase at 1000 K and l.'bgl cm? ^ 4 data points from Ice 
X phase at 2000 K and 'i.2'bg/av? to A.QQg / civ? ^ and the rest of the data points from the 
superionic phase at 4000 K and ^.QQg/civ? up to Ibg/cm?. Since the temperature effect on 
density becomes smaller towards higher pressure, all the isothermal pressure- density curves 
converge on the isentropic pressure-density curves as well as the pressure-density curve along 
the melting curve. 

The FMNR2009's EOS has been confirmed experimentally by Thomas Mattson et al. 
at the Sandia National Laboratories. At 8.893 TPa, FMNR2009's EOS is switched to the 



TED EOS in lSeaeer. Kuchner. Hier-Maiumder. k Militzed 120071 (SKHMM2007). 



3.34. SKHMM2007's EOS 



At ultrahigh pressure, the effect of electron-electron interaction can be safely ignored and 
electrons can be treated as a gas of non-interacting particles that obey the Pauli exclusion 
principle subject to the Coulomb field of the nuclei. Assuming the Coulomb potential is 
spatially slowly varying throughout the electron gas that the electronic wave functions can 
be approximated locally as plane waves, the so-called TED so lution could be derived so that 
the P auli exclusion pressure balances out the Coulomb forces (lEliezer et al.ll2002[ iMacfarlane 

igsJ). 



In SKHMM2007, a modified TED bv lSalpeter fc Zapolskvl fll967f ) is used. It is modified 
in the sense that the authors have added in a density- dependent correlation energy term 
which characterizes electron interaction effects. 



Here, ISeaeer et all 's EOS is slightly modified to connect to the EMNR2 009's EOS. A t 



8.893 TPa, the ratio of the density p between the EMN R2009's EOS and ISeager et al. 



EOS is 1.04464. This ratio is multiplied to the original ISeager et al.r s EOS density p to 
produce the actual EOS used in our calculation for p> 8.893 TPa. Only the most massive 
planets in our model could reach this pressure in the Il20-layer, so this choice of EOS for 
p> 8.893 TPa has small effect on the overall mass-radius relation to be discussed in the next 
section. 



The EOS of H2O is shown in Fig. [T]as the lower curve (blue curve). 
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4. Result 

4.1. Mass-Radius Contours 

Given mass and radius input, various sets of mass-radius contours can be used to quickly 
determine the characteristic interior structure quantities of a 2-layer planet including its pO 
(central pressure), pl/pO (ratio of core-mantle boundary pressure over central pressure), 
CMF (core mass fraction), and CRF (core radius fraction). 

The 2-layer model is uniquely solved and represented as a point on the mass-radius 
diagram given any pair of two parameters from the following list: M (mass), R (radius), pO, 
pl/pO, CMF, CRF, etc. The contours of constant M or R are trivial on the mass-radius 
diagram, which are simply vertical or horizontal lines. The contours of constant pO, pl/pO, 
CMF, or CRF are more useful. 

Within a pair of parameters, fixing one and continuously varying the other, the point on 
the mass-radius diagram moves to form a curve. Multiple values of the fixed parameter give 
multiple parallel curves forming a set of contours. The other set of contours can be obtained 
by exchanging the fixed parameter for the varying parameter. The two sets of contours 
crisscross each other to form a mesh, which is a natural coordinate system (see Fig. |2]) of 
this pair of parameters, superimposed onto the existing Cartesian (M, R) coordinates of the 
mass-radius diagram. This mesh can be used to determine the two parameters given the 
mass and radius input or vice versa. 

4.1.1. Fe-MgSiOs planet 

As an example, the mesh of pO with pl/pO for Fe-core MgSiOa-mantle planet is illus- 
trated in the first subplot (upper left corner) of Fig. |2l The mesh is formed by pO-contours 
and pl/pO-contours crisscrossing each other. The more vertical set of curves represents the 
pO-contours. The ratio of adjacent pO-contours is 10°'^ (see Tabled]). The more horizontal 
set of curves represents the pl/pO-contours. From bottom up, the pl/pO values vary from 
to 1 with step size 0.1. 

Given a pair of pO and pi input, users may interpolate from the mesh to find the 
mass and radius. On the other hand, given the mass and radius of a planet, users may 
also interpolate from the mesh to find the corresponding pO and pi of a 2-layer Fe core 
MgSiOa-mantle planet. 

Similarly, the contour mesh of pO with CMF for the Fe-MgSiOa planet is shown as the 
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Fig. 2. — Mass-Radius contours of 2-layer planet. 1st row: Fe-MgSiOa planet. 2nd row: 
MgSi03-H20 planet. 3rd row: Fe-H20 planet. 1st column: contour mesh of pl/pO with pO. 
2nd column: contour mesh of CMF with pO. 3rd column: contour mesh of CRF with pO. To 
find out what pO value each pO-contour corresponds to, please refer to Table [H 
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second subplot from the left in the first row of Fig. [2l As a reference point, for a pure-Fe 
planet with pO = 10^^ Pa, M = 0.1254Me, R = 0.417/?®. 

The contour mesh of pO with CRF for the Fe-MgSiOs planet is shown as the third 
subplot from the left of the first row of Fig. |2j 



4.1.2. MgSiOi-H20 planet 

For 2-layer MgSi03-H20 planet, the 3 diagrams (pO contours pair with pl/pO contours, 
CMF contours, or CRF contours) are the subplots of the second row of Fig. [2J 

As a reference point, for a pure-MgSiOs planet with = 10^°'^ Pa, M = O.122M0, R = 
0.5396i?ffi. 



4.1.3. Fe-H2 planet 

For 2-layer Fe-H20 planet, the 3 diagrams (pO contours pair with pl/pO contours, CMF 
contours, or CRF contours) are the subplots of the third row of Fig. |2l 

As a reference point, for a pure-Fe planet with pO = 10^^ Pa, M = 0.1254:M^, R = 
0.417i?ffi. 



4.2. Mass-Radius Curves 



For observers' interest, 6 characteristic mass-radius curves are plotted (Fig. |3]) and tab- 
ulated (Tabled]), representing the pure-Fe planet, half-Fe half-MgSiOs planet, pure MgSiOs 
planet, half-MgSiOa half-HaO planet, 75% H20-25% MgSiOg planet, and pure H2O planet. 
These fractions are mass fractions. Fig. [3] also shows some recently discovered exoplan- 
ets within the relevant mass-radius r egime for compariso n. These pl anets include Kepler - 
10b JSatalha et al.lboill). K epler-llb Jhissauer et al.lboill) Kepler-llf d hissauer et al"l201lh. 
Kepler-18b dCochran et aDboilh . Kepler-36b dCarter et al.ll2012[). and KeDl er-20b.c.d dlll et al, 
2012h . They also include Kepler-20e {R = 0.868+g:gI^i?e dPressin et"ai] 120121 ) . the mass 



range i s determined by pur e-silicate mass-radius curve and the maximum coUision al stripping 
curve dMarcus et allboioh ). Kepler-20f {R = 1.034+°;}°°/?® JPressin et all 120121 ) . the mass 
range is determined by 75 % water-ice and 25% silicate mass-radius curves and t he maximum 
coUis ional stripping curve ( iMarcus et al. II2OI0I)), Kepler-21b (R = 1.64±0.04i?e dHowell et all 
2OI2I ). The upper limit for mass is 10. 4M®: the 2-cr upper limit preferred in the paper. The 



- 14 - 



lower limit is 4M®, which is in between the "Earth" and "50% H2O-50% MgSiOs" model 
curves - the planet is very hot and is unlikely to h ave much water content if any at all.), 
Kepler-22b {R = 2.38 ±0.13i?e JSorucki et al.ll2012h . The 1-a upper limit for mass is 36M® 
for an eccentri c orbit (or 27Mm for circular orbit)) CoRoT-7b (M = 7.42 ± 1.21Me, R = 



2009 



Queloz et al. 



20091)) 55 Cancri e (M 



1.58 ± O.li?© flHatzes et al.ll2011l : iLeeer et al 

8.63 ± 0.35Me, R = 2.00 ± O.UR^ (1win~ al.ll201lh ). and GJ 1214b fICharbonneau et al. 



20091). 



4.3. Levels of planet differentiation: the effect of Fe partitioning between 

mantle and core 

All models of planets discussed so far assume that all Fe is in the core, while all Mg, 
Si and O are in the mantle, i.e., that a planet is fully differentiated. However, we know 
that in terrestrial planets some Fe is incorporated into the mantle. There are two separate 
processes which affect the Fe content of the mantle: (1) mechanical segregation of Fe-rich 
metal from the mantle to the core, and (2) different redox conditions resulting in a different 
Fe/Mg ratio within the mantle, which in turn affects the relative size of the core and mantle. 
In this section we show the effects of (1) undifferentiated versus fully-differentiated, and (2) 
reduced versus oxidized planetary structure on the mass-radius relation for a planet with the 
same Fe/Si and Mg/Si ratios. 

For simplicity, here we ignore the H2O and gaseous content of the planet and only 
consider the planet made of Fe, Mg, Si, O. To facilitate comparison between different cases, 
we fix the global atomic ratios of Fe/Mg=l a nd Mg/Si=l; these fit well within the range of 



local stellar abundances ( iGrasset et al.l 120091 ). 



In particular, we consider a double-layer planet with a core and a mantle in two sce- 
narios. One follows the incomplete mechanical separation of the Fe-rich metal during planet 
formation, and results in addition of Fe to the mantle as metal particles. It does not change 
EOS of the silicate components, but requires adding an Fe-EOS to the mantle mixture. Thus, 



^mass in Earth Mass (Af® = 5.9742 x lO^^fcg) 
^radius in Earth Radius (i?® — 6.371 x lO^m) 

■^CRF stands for core radius fraction, the ratio of the radius of the core over the total radius of the planet, 
in the two-layer model 

''pl/pO stands for core-mantle-boundary pressure fraction, the ratio of the pressure at the core-mantle 
boundary (pi) over the pressure at the center of the planet (pO), in the two-layer model 
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5 10 

M (Earth Mass) 



Fig. 3. — Currently known transiting exoplanets are shown with their measured mass and 
radius with observation uncertainties. Earth and Venus are shown for comparison. The 
curves are calculated for planets composed of pure Fe, 50% Fe-50% MgSiOs, pure MgSiOs, 
50% H2O-50% MgSiOs, 75% H20-25% MgSiOg an d pure H^O. The red dashed curve is the 
maximum collisional stripping curve calculated by iMarcus et al.l ( 120101 ) . 



Table 1: Data for the 6 characteristic mass-radius curves 



50%Fe-50%MgSiO:i 



MgSiOs 



50%MgSiO3-50%H2O 



25%MgStO:r75%H20 



H2O 



logn){p{)) 


Mass' 


Radius^ 


Mass 


Radius 


9.6 


0.001496 


0.09947 


0.00177 


0.121 


9.7 


0.002096 


0.1112 


0.002481 


0.1354 


9.8 


0.002931 


0.1243 


0.003472 


0.1513 


9.9 


0.00409 


0.1387 


0.004848 


0.169 


10 


0.005694 


0.1546 


0.006752 


0.1885 


10.1 


0.007904 


0.1722 


0.00938 


0.2101 


10.2 


0.01094 


0.1915 


0.01299 


0.2339 


10.3 


0.01507 


0.2126 


0.01792 


0.26 


10.4 


0.0207 


0.2356 


0.02464 


0.2886 


10.5 


0.02829 


0.2606 


0.03372 


0.3197 


10.6 


0.0385 


0.2876 


0.04595 


0.3535 


10.7 


0.05212 


0.3167 


0.06231 


0.3901 


10.8 


0.07021 


0.348 


0.08408 


0.4296 


10.9 


0.09408 


0.3814 


0.1129 


0.4721 


11 


0.1254 


0.417 


0.1508 


0.5177 


11.1 


0.1663 


0.4548 


0.2003 


0.5663 


11.2 


0.2193 


0.4949 


0.2647 


0.618 


11.3 


0.2877 


0.537 


0.348 


0.6728 


11.4 


0.3754 


0.5814 


0.455 


0.7307 


11.5 


0.4875 


0.6279 


0.5919 


0.7916 


11.6 


0.6298 


0.6765 


0.7659 


0.8554 


11.7 


0.8096 


0.727 


0.986 


0.9221 


11.8 


1.036 


0.7796 


1.261 


0.9907 


11.9 


1.319 


0.834 


1.606 


1.062 


12 


1.671 


0.8902 


2.036 


1.136 


12.1 


2.108 


0.9481 


2.568 


1.211 


12.2 


2.618 


1.007 


3.226 


1.289 


12.3 


3.31 


1.068 


4.032 


1.369 


12.4 


4.119 


1.13 


5.018 


1.451 


12.5 


5.103 


1.193 


6.216 


1.534 


12.6 


6.293 


1.257 


7.665 


1.618 


12.7 


7.727 


1.321 


9.39 


1.7 


12.8 


9.445 


1.386 


11.45 


1.783 


12.9 


11.49 


1.451 


13.91 


1.865 


13 


13.92 


1.515 


16.81 


1.947 


13.1 


16.78 


1.579 


20.21 


2.027 


13.2 


20.14 


1.642 


24.2 


2.106 


13.3 


24.05 


1.704 


28.85 


2.183 


13.4 


28.6 


1.765 


34.23 


2.258 


13.5 


33.87 


1.824 


40.45 


2.331 


13.6 


39.94 


1.881 


47.59 


2.401 


13.7 


46.92 


1.937 


55.76 


2.468 


13.8 


54.93 


1.99 


65.07 


2.531 


13.9 


64.08 


2.042 


75.65 


2.59 


14 


74.51 


2.091 


87.65 


2.645 


14.1 


86.37 


2.138 


101.2 


2.697 


14.2 


99.8 


2.183 


116.5 


2.745 


14.3 


115 


2.226 


133.8 


2.789 


14.4 


132.1 


2.266 


153.1 


2.829 



CRF 



pl/pO' 


Mass 


Radius 


Mass 


Radius 


CRF 


0.2931 


0.00623 


0.2029 


0.008278 


0.2963 


0.5969 


0.2925 


0.008748 


0.227 


0.01156 


0.3286 


0.6013 


0.2918 


0.01227 


0.254 


0.01615 


0.3647 


0.6051 


0.2909 


0.01717 


0.2838 


0.02255 


0.4049 


0.6085 


0.2898 


0.02399 


0.317 


0.0313 


0.4484 


0.6119 


0.2886 


0.03343 


0.3536 


0.04314 


0.4944 


0.6165 


0.2871 


0.04644 


0.3939 


0.05943 


0.5449 


0.6209 


0.2855 


0.0643 


0.4381 


0.0817 


0.6003 


0.6249 


0.2836 


0.08866 


0.4865 


0.112 


0.6607 


0.6285 


0.2815 


0.1217 


0.5391 


0.1528 


0.7262 


0.6318 


0.2792 


0.1661 


0.596 


0.2074 


0.7966 


0.6347 


0.2767 


0.2255 


0.6573 


0.2799 


0.8718 


0.6373 


0.274 


0.3042 


0.7228 


0.3754 


0.9515 


0.639C 


0.2711 


0.4075 


0.7925 


0.4998 


1.035 


0.6416 


0.2682 


0.542 


0.8661 


0.6607 


1.123 


0.6434 


0.2651 


0.7143 


0.9429 


0.8653 


1.214 


0.6448 


0.2621 


0.927 


1.02 


1.117 


1.305 


0.6455 


0.259 


1.2 


1.101 


1.437 


1.4 


0.6458 


0.2559 


1.545 


1.186 


1.842 


1.499 


0.6458 


0.2529 


1.981 


1.274 


2.351 


1.602 


0.6453 


0.25 


2.525 


1.365 


2.988 


1.708 


0.6445 


0.2473 


3.203 


1.458 


3.78 


1.818 


0.6435 


0.2451 


4.043 


1.554 


4.763 


1.933 


0.6423 


0.2429 


5.077 


1.653 


5.972 


2.05 


0.641 


0.2407 


6.297 


1.749 


7.392 


2.165 


0.6394 


0.2385 


7.714 


1.842 


9.043 


2.275 


0.6372 


0.2363 


9.423 


1.935 


11.03 


2.386 


0.6345 


0.2342 


11.47 


2.029 


13.4 


2.498 


0.6317 


0.2321 


13.87 


2.121 


16.18 


2.608 


0.6285 


0.23 


16.73 


2.213 


19.48 


2.717 


0.6252 


0.228 


20.1 


2.304 


23.36 


2.825 


0.6219 


0.2267 


24.07 


2.394 


27.94 


2.933 


0.619 


0.2254 


28.68 


2.483 


33.24 


3.037 


0.6162 


0.2241 


33.99 


2.569 


39.33 


3.138 


0.6134 


0.2227 


40.05 


2.65 


46.26 


3.234 


0.6105 


0.2213 


46.88 


2.728 


54.07 


3.325 


0.6075 


0.2198 


54.49 


2.799 


62.77 


3.409 


0.6042 


0.2182 


63.08 


2.866 


72.58 


3.488 


0.6006 


0.2166 


72.75 


2.928 


83.62 


3.563 


0.5967 


0.2151 


83.59 


2.986 


95.97 


3.631 


0.5928 


0.2136 


95.68 


3.039 


109.8 


3.695 


0.5888 


0.2123 


109.1 


3.088 


125.1 


3.754 


0.5847 


0.2114 


124 


3.132 


142 


3.807 


0.5807 


0.2106 


140.3 


3.17 


160.6 


3.854 


0.5768 


0.21 


158.2 


3.204 


181 


3.895 


0.5728 


0.2094 


177.8 


3.232 


203.3 


3.93 


0..5688 


0.2089 


199.2 


3.255 


227.7 


3.959 


0.5648 


0.2085 


222.5 


3.274 


254.2 


3.983 


0..5fi07 


0.2082 


248.1 


3.288 


283.3 


4.002 


0.5567 



pl/pO 



Mass 


Radius 


CRF 


pl/pO 


Mass 


Radius 


0.01217 


0.3616 


0.4413 


0.3782 


0.04663 


0.58 


0.01699 


0.4009 


0.4444 


0.3803 


0.06174 


0.63 


0.02351 


0.4432 


0.4475 


0.385 


0.08208 


0.6853 


0.03213 


0.4869 


0.4513 


0.3944 


0.1091 


0.7455 


0.04399 


0.535 


0.4554 


0.4016 


0.1445 


0.8102 


0.06029 


0.588 


0.4592 


0.4073 


0.1904 


0.879 


0.08255 


0.6461 


0.4629 


0.4119 


0.2494 


0.9515 


0.1127 


0.7093 


0.4662 


0.4159 


0.3249 


1.028 


0.1533 


0.7777 


0.4693 


0.4197 


0.4206 


1.107 


0.2074 


0.8509 


0.4721 


0.4234 


0.5416 


1.19 


0.2789 


0.9289 


0.4747 


0.427 


0.6938 


1.276 


0.3726 


1.011 


0.477 


0.4305 


0.8866 


1.366 


0.4946 


1.098 


0.4792 


0.4339 


1.132 


1.461 


0.6517 


1.188 


0.4811 


0.437 


1.444 


1.562 


0.8529 


1.282 


0.4827 


0.4397 


1.841 


1.669 


1.107 


1.38 


0.4839 


0.4411 


2.346 


1.782 


1.42 


1.477 


0.4843 


0.44 


2.985 


1.901 


1.818 


1.58 


0.484 


0.44 


3.77 


2.023 


2.321 


1.688 


0.4834 


0.4403 


4.735 


2.147 


2.957 


1.802 


0.4828 


0.4399 


5.909 


2.274 


3.752 


1.92 


0.4822 


0.4391 


7.325 


2.401 


4.732 


2.041 


0.4813 


0.4383 


9.038 


2.529 


5.936 


2.164 


0.4803 


0.4377 


11.11 


2.66 


7.407 


2.289 


0.4792 


0.4373 


13.55 


2.789 


9.127 


2.411 


0.4778 


0.4341 


16.42 


2.915 


11.12 


2.529 


0.4755 


0.4299 


19.77 


3.039 


13.52 


2.647 


0.4727 


0.4275 


23.68 


3.16 


16.37 


2.765 


0.4695 


0.4265 


28.21 


3.278 


19.72 


2.882 


0.467 


0.4248 


33.49 


3.393 


23.68 


2.998 


0.4646 


0.4237 


39.62 


3.506 


28.31 


3.111 


0.4623 


0.423 


46.72 


3.616 


33.71 


3.222 


0.46 


0.4225 


54.92 


3.724 


39.97 


3.33 


0.4579 


0.4214 


64.22 


3.826 


47.15 


3.434 


0.4556 


0.4199 


74.79 


3.924 


55.31 


3.534 


0.4531 


0.4177 


86.85 


4.017 


64.47 


3.627 


0.4503 


0.4148 


100.3 


4.104 


74.65 


3.713 


0.4472 


0.4114 


115.3 


4.183 


86.14 


3.793 


0.4438 


0.4087 


131.9 


4.256 


99.08 


3.87 


0.4404 


0.4064 


150.3 


4.322 


113.6 


3.94 


0.437 


0.4043 


170.8 


4.382 


129.7 


4.006 


0.4338 


0.4022 


193.6 


4.435 


147.6 


4.064 


0.4306 


0.4001 


218.7 


4.483 


167.2 


4.116 


0.4274 


0.3979 


246.6 


4.525 


188.8 


4.162 


0.4242 


0.3956 


277.3 


4.561 


212.5 


4.202 


0.4209 


0.3932 


311.3 


4.592 


238.4 


4.236 


0.4175 


0.3908 


348.7 


4.618 


266.8 


4.265 


0.4142 


0.3884 


390.1 


4.639 


297.9 


4.288 


0.4109 


0.3861 


435.9 


4.656 


332.1 


4.307 


0.4076 


0.3841 


486.4 


4.669 



0.6892 
0.6888 
0.6882 
0.6876 
0.6868 
0.6858 
0.6847 
0.6834 
0.6819 
0.6801 
0.6782 
0.6759 
0.6735 
0.6708 
0.6679 
0.6648 
0.6616 
0.6582 
0.6547 
0.6512 
0.6476 

0.644 
0.6407 
0.6374 
0.6342 

0.631 
0.6278 
0.6247 
0.6216 
0.6184 
0.6153 
0.6126 

0.61 
0.6074 
0.6048 
0.6024 
0.5999 
0.5973 
0.5947 
0.5921 
0.5895 
0.5869 
0.5844 

0.582 
0.5797 
0.5775 
0.5755 
0.5735 
0.5716 



0.2374 
0.2398 
0.2413 
0.242 
0.2446 
0.2495 
0.253 
0.2554 
0.2571 
0.2585 
0.2595 
0.2603 
0.261 
0.2615 
0.2618 
0.2615 
0.2598 
0.2589 
0.2582 
0.2574 
0.2565 
0.2552 
0.2536 
0.252 
0.2488 
0.2449 
0.242 
0.2396 
0.2371 
0.2353 
0.2339 
0.2323 
0.2306 
0.2288 
0.2266 
0.2238 
0.2207 
0.2178 
0.2152 
0.2129 
0.2106 
0.2083 
0.2061 
0.2038 
0.2015 
0.1992 
0.1971 
0.1951 
0.1932 
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the planet generally consists of a Fe metal core and a partially differentiated mantle consist- 
ing of the mixture of metallic Fe and MgSiOs silicates. While the distribution of metallic 
Fe may have a radial gradient, for simplicity we assume that it is uniformly distributed in 
the silicate mantle. Within scenario 1, we calculate three cases to represent different levels 
of differentiation: 



Case 1: complete differentiation metallic Fe core and MgSiOa silicate mantle. ForFe/Mg=l, 
CMF=0.3574. 

Case 2: partial differentiation half the Fe forms a smaller metallic Fe core, with the 
other half of metal being mixed with MgSiOs silicates in the mantle. CMF=0.1787. 

Case 3: no differentiation All the metallic Fe is mixed with MgSiOs in the mantle. 
CMF=0 fno core). 



The other scenario assumes different redox conditions, resulting in different Fe/Mg ratios 
in mantle minerals, and therefore requiring different EOS for the Fe^"'"-bearing silicates and 
oxides. More oxidized mantle means adding more Fe in the form of FeO to the MgSiOs 
silicates to form (Mg,Fe)Si03 silicates and (Mg,Fe)0 magnesiowiistite (mv), thus reducing 
the amount of Fe in the core. The exact amounts of (Mg,Fe)Si03 and (Mg,Fe)0 in the 
mantle are determined by the following mass balance equation: 



xFeO + MgSi03 ^ {Mg Fe^)Si03 + x {Mg^, Fe^)0 (6) 

X denotes the relative amount of FeO added to the silicate mantle, x=0 being the most 
reduced state with no Fe in the mantle, and x=l being the most oxidized state with all Fe 
existing as oxides in the mantle. This oxidization process conserves the global Fe/Mg and 
Mg/Si ratios, but increases O content and thus the 0/Si ratio of the planet since Fe is added 
to the mantle in the form of FeO. Because in stellar environments O is excessively abun dant 



relative to Mg, Si, and Fe (e.g., the solar elemental abundances ( lAsplund et al.ll2009l )). it 
is not a limiting factor in our models of oxidized planets. We calculate the following three 
cases to represent the full range of redox conditions: 



Case 4: no oxidization of Fe x=0. Metal Fe core and MgSiOa silicate mantle. For 
Fe/Mg=l, 0/Si=3, CMF=0.3574. 

Case 5: partial oxidization of Fe x=0.5. Half the Fe forms smaller metal core, the other 
half is added as FeO to the mantle. 0/Si=3.5, CMF=0.1700. 
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Case 6: complete oxidization of Fe x=l. All Fe is added as FeO to the mantle, result- 
ing in no metal core at all. 0/Si=4, CMF=0. 



Notice that Case 4 looks identical to Case 1. However, the silicate EOS used to calcu- 
late Case 4 is different from Case 1 at ultrahigh pressures (beyond .90 T Pa). For Cases 4, 
5, and 6, the (Mg,Fe)Si03 EOS is adopted from [Caracas fc CohenI (120081 ) and ISpera et al. 



( 120061 ) which only consider perovskite (pv) and post-perovskite (ppv) phases without includ- 
ing further dissociation beyond 0.90 TPa, since the Fe^"'"-bearing silicate EOS at ultrahigh 
pressures is hardly available. On the other hand. Comparison between Case 1 and Case 4 
also shows the uncertainty on mass-radius relation resulting from the different choice of EOS 
(see Table E]). Fe^^-bearing pv and ppv have the general formula: (Mgy,Fei_y)Si03, where 
y denotes the relative atomic number fraction of Mg and Fe in the silicate mineral. The 
(Mgy,Fei_j^)Si03 sihcate is therefore a binary component equilibrium solid s olution. It could 
either be pv or ppv or both depending on the pressure (ISpera et al.ll2006l ). We can safely 
approximate the narrow pressure region where pv and ppv co-exist as a single transition 
pressure from pv to ppv. This pressure is calculated as the arithmetic mean of the initial 
transition pr essure of pv — )■ pv+ ppv mixture and the final transition pressure of pv+ppv mix- 
ture —7- ppv (ISpera et al.l 120061 ). The pv EOS and ppv EOS are connected at this transition 
pressure to form a consistent EOS for all pressures. 

Addition of FeO to MgSiOa results in the appearance of a second phase, magnesiowiistite, 
in th e mantle accord ing to Eq. |6l The (Mg,Fe)0 EOS for Cases 4, 5, and 6 is adopted 
from iFei et al.l (120071 ). which includes the electronic spin transition of high-spin to low-spin 
in Fe^"*". For simplicity, we assume that pv/ppv and mw have the same Mg/Fe ratio. 

Fig. m shows fractional differences [rj) in radius of Cases 2 & 3 compared to Case 1 as 
well as Cases 5 & 6 compared to Case 4 (ro is radius of the reference case, which is that of 
Case 1 for Cases 2 & 3 and is that of Case 4 for Case 5 & 6): 



V 



r — vq 
ro 



(7) 



Oxidization of Fe (partitioning Fe as Fe-oxides from the core into the mantle) makes the 
planet appear larger. The complete oxidization of Fe makes the radius 3% larger for small 
planets around 1 M®, then the difference decreases with increasing mass within the mass 
range of 1 to 20 M^. Undifferentiated planets (partitioning of metallic Fe from the core into 
the mantle) appear smaller than fully differentiated planets. The completely undifferentiated 
planet is practically indistinguishable in radius for small planets around 1 M®, then the 
difference increases to 1%-level around 20 M®. The mass, radius, CRF and pl/pO data of 
Cases 1 through 6 are listed in Table [2J 
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Case 6: complete oxidization of Fe 



Case 5: partial (50%) oxidization of Fe 



Case 1 & Case 4: differentiated & no oxidization of Fe 



Case 2: partially (50%) differentiated 




5 10 15 

M (Earth Mass) 



20 



Fig. 4. — fractional differences in radius resulting from Fe partitioning between mantle and 
core, black curve: Case 1 & Case 4 (complete differentiation and no oxidization of Fe); solid 
brown curve: Case 2 (partial (50%) differentiation: 50% metallic Fe mixed with the mantle); 
solid pink curve: Case 3 (no differentiation: all metallic Fe mixed with the mantle); dashed 
brown curve: Case 5 (partial (50%) oxidization of Fe); dashed pink curve: Case 6 (complete 
oxidization of Fe) 



Table 2: Data of Cases 



1 through 6 



Case 1 (CMF=0.3574) 



Case 2 (CMF=0.1787) 



Case 3 



Case 4 (CMF=0.3574) 



Case 5 (CMF=0.1700) 



Case 6 



Mass' Radius^ CRF'' pl/pO* 



Mass Radius CRF pl/pO 



Ma 



Radius 



Mass Radius CRF 



pl/pO 


Mass 


Radius 


CRF 


pl/pO 


Mass 


Radius 


0.3847 


0.002526 


0.1428 


0.4584 


0.5625 


0.004921 


0.1804 


0.3839 


0.003542 


0.1597 


0.4579 


0.5617 


0.006905 


0.2018 


0.383 


0.004959 


0.1786 


0.4574 


0.5607 


0.009673 


0.2256 


0.382 


0.00693 


0.1995 


0.4567 


0.5596 


0.01352 


0.252 


0.3807 


0.009662 


0.2227 


0.456 


0.5582 


0.01887 


0.2812 


0.3792 


0.01344 


0.2484 


0.455 


0.5567 


0.02625 


0.3135 


0.3774 


0.01863 


0.2766 


0.454 


0.5548 


0.03639 


0.3489 


0.3754 


0.02574 


0.3076 


0.4527 


0.5528 


0.05027 


0.3878 


0.3732 


0.03544 


0.3416 


0.4513 


0.5504 


0.06915 


0.4301 


0.3706 


0.04859 


0.3786 


0.4497 


0.5479 


0.09465 


0.476 


0.3679 


0.06631 


0.4189 


0.448 


0.5451 


0.1289 


0.5257 


0.3649 


0.09006 


0.4624 


0.446 


0..5422 


0.1744 


0.5789 


0.3617 


0.1217 


0.5094 


0.4439 


0.5392 


0.2345 


0.6358 


0.3584 


0.1636 


0.5599 


0.4417 


0.5361 


0.3113 


0.6946 


0.3-55 


0.2187 


0.6139 


0.4394 


0.533 


0.408 


0.7-548 


0.3515 


0.2903 


0.6709 


0.4371 


0.5302 


0..5315 


0.8177 


0.348 


0.3833 


0.7313 


0.4347 


0.5275 


0.6905 


0.884 


0.3445 


0.5024 


0.7945 


0.4323 


0.5255 


0.8932 


0.9536 


0.3412 


0.6542 


0.8604 


0.43 


0.5237 


1.15 


1.026 


0.3379 


0.8478 


0.9294 


0.4278 


0.5219 


1.472 


1.101 


0.3352 


1.093 


1.001 


0.4255 


0.5201 


1.874 


1.178 


0.333 


1.402 


1.076 


0.4233 


0.5184 


2.373 


1.258 


0.3307 


1.789 


1.154 


0.4211 


0.5167 


2.99 


1.34 


0.3284 


2.272 


1.233 


0.419 


0.5152 


3.75 


1.423 


0.3262 


2.871 


1.316 


0.4168 


0.5137 


4.68 


1.509 


0.324 


3.61 


1.4 


0.4147 


0.5123 


5.818 


1.596 


0.3219 


4.519 


1.486 


0.4126 


0.5109 


7.203 


1.684 


0.3197 


5.63 


1.574 


0.4105 


0.5095 


8.886 


1.775 


0.3176 


6.983 


1.664 


0.4085 


0.5081 


10.93 


1.866 


0.3155 


8.626 


1.754 


0.4064 


0..5065 


13.39 


1.96 


0.3133 


10.01 


1.847 


0.4043 


0.5049 


16.37 


2.054 


0.311 


13 


1.94 


0.4021 


0.503 


19.96 


2.15 


0.3086 


15.86 


2.034 


0.3999 


0..5009 


24.29 


2.248 


0.3061 


19.29 


2.129 


0.3976 


0.4985 


29.48 


2.347 


0.3033 


23.36 


2.224 


0.3952 


0.4958 


35.73 


2.448 


0.3004 


28.19 


2.32 


0.3927 


0.4928 


43.22 


2.551 


0.2972 


33.9 


2.417 


0.3901 


0.4893 


52.21 


2.656 


0.2937 


40.62 


2.513 


0.3873 


0.4853 


62.98 


2.762 


0.29 


48.52 


2.61 


0.3844 


0.481 


75.88 


2.87 


0.2862 


57.81 


2.708 


0.3813 


0.4765 


91.32 


2.981 


0.2821 


68.68 


2.806 


0.3782 


0.4714 


109.8 


3.093 


0.2779 


81.38 


2.905 


0.375 


0.466 


131.9 


3.208 


0.2735 


96.2 


3.004 


0.3716 


0.4604 


158.4 


3.326 


0.2689 


113.5 


3.104 


0.3682 


0.4544 


190 


3.446 


0.2642 


133.6 


3.204 


0.3647 


0.4481 


227.8 


3.568 


0.2593 


157 


3.306 


0.3612 


0.4415 


273 


3.694 


0.2542 


184.1 


3.408 


0.3575 


0.4345 


327.1 


3.823 


0.2489 


215.5 


3.511 


0.3538 


0.4271 


391.8 


3.955 


0.2435 


251.8 


3.615 


0.3499 


0.4195 


469.1 


4.091 



9.6 


0.002009 


0.1302 


9.7 


0.002816 


0.1456 


9.8 


0.003941 


0.1628 


9.9 


0.005504 


0.1819 


10 


0.00767 


0.203 


10.1 


0.01066 


0.2263 


10.2 


0.01477 


0.252 


10.3 


0.02039 


0.2802 


10.4 


0.02805 


0.311 


10.5 


0.03842 


0.3447 


10.6 


0.05239 


0.3814 


10.7 


0.07111 


0.4211 


10.8 


0.09605 


0.464 


10.9 


0.1291 


0.5101 


11 


0.1725 


0.5595 


11.1 


0.2294 


0.6123 


11.2 


0.3033 


0.6684 


11.3 


0.399 


0.7278 


11.4 


0.5219 


0.7906 


11.5 


0.679 


0.8565 


11.6 


0.8779 


0.9251 


11.7 


1.128 


0.996 


11.8 


1.443 


1.07 


11.9 


1.837 


1.146 


12 


2.327 


1.225 


12.1 


2.934 


1.306 


12.2 


3.682 


1.39 


12.3 


4.6 


1.475 


12.4 


5.72 


1.562 


12.5 


7.071 


1.649 


12.6 


8.686 


1.736 


12.7 


10.62 


1.822 


12.8 


12.93 


1.908 


12.9 


15.67 


1.993 


13 


18.9 


2.078 


13.1 


22.69 


2.161 


13.2 


27.12 


2.242 


13.3 


32.27 


2.322 


13.4 


38.2 


2.399 


13.5 


45.05 


2.472 


13.6 


52.85 


2.542 


13.7 


61.72 


2.608 


13.8 


71.81 


2.67 


13.9 


83.26 


2.728 


14 


96.2 


2.781 


14.1 


110.8 


2.831 


14.2 


127.2 


2.876 


14.3 


145.5 


2.917 


14.4 


166 


2.954 



0.5972 
0.5967 
0.5961 
0.5954 
0.5946 
0.5936 
0.5924 
0.591 
0.5894 
0.5876 
0.5856 
0.5833 
0.5809 
0.5782 
0.5753 
0.5723 
0.5692 
0.566 
0.5627 
0.5594 
0.5561 
0.553 

0.55 
0.547 
0.5441 
0.5412 
0.5383 
0.5355 
0.5327 

0.-53 
0.5275 
0.5251 
0.5228 
0.5206 
0.5183 
0.5161 
0.5139 
0.5116 
0.5093 
0.5069 
0.5046 
0.5024 
0.5002 
0.4982 
0.4962 
0.4944 
0.4926 
0.491 
0.4894 



0.3847 
0.3839 

0.383 

0.382 
0.3807 
0.3792 
0.3774 
0.3754 
0.3732 
0.3706 
0.3679 
0.3649 
0.3617 
0.3584 

0.355 
0.3515 

0.348 
0.3445 
0.3412 
0.3379 
0.3352 

0.333 
0.3307 
0.3284 
0.3262 

0.324 
0.3219 
0.3197 
0.3176 
0.3161 
0.3152 
0.3143 
0.3134 
0.3125 
0.3116 
0.3104 
0.3092 
0.3078 
0.3066 
0.3056 
0.3048 
0.3044 
0.3042 
0.3041 
0.3042 
0.3043 
0.3045 
0.3049 
0.3054 



0.002395 


0.1381 


0.4736 


0.003359 


0.1545 


0.4731 


0.004704 


0.1727 


0.4726 


0.006573 


0.1929 


0.4719 


0.009165 


0.2154 


0.4711 


0.01275 


0.2401 


0.4702 


0.01767 


0.2674 


0.469 


0.02442 


0.2974 


0.4678 


0.03362 


0.3303 


0.4663 


0.0461 


0.3661 


0.4647 


0.06292 


0.405 


0.4628 


0.08548 


0.4472 


0.4608 


0.1155 


0.4927 


0.4.587 


0.15.53 


0.5416 


0.4564 


0.2077 


0.5939 


0.454 


0.2762 


0.6495 


0.4515 


0.3652 


0.7085 


0.4489 


0.48 


0.7708 


0.4463 


0.6266 


0.8358 


0.4437 


0.8118 


0.9031 


0.4412 


1.047 


0.9733 


0.4388 


1.343 


1.046 


0.4365 


1.715 


1.122 


0.4341 


2.179 


1.2 


0.4318 


2.755 


1.281 


0.4295 


3.468 


1.364 


0.4273 


4.344 


1.448 


0.425 


5.403 


1.533 


0.4229 


6.671 


1.617 


0.4208 


8.202 


1.702 


0.4189 


10.04 


1.787 


0.417 


12.23 


1.871 


0.4152 


14.83 


1.955 


0.4134 


17.9 


2.038 


0.4116 


21.53 


2.12 


0.4098 


25.76 


2.201 


0.4079 


30.67 


2.279 


0.406 


36.33 


2.354 


0.404 


42.8 


2.426 


0.402 


50.18 


2.494 


0.4001 


58.56 


2.557 


0.3982 


68.09 


2.618 


0.3963 


78.88 


2.674 


0.3945 


91.07 


2.727 


0.3928 


104.8 


2.775 


0.3911 


120.1 


2.819 


0.3895 


137.3 


2.859 


0.388 


156.4 


2.894 


0.3865 


177.5 


2.925 


0.3851 



0.5634 
0.5626 
0.5616 
0.5604 
0.5591 
0.5575 
0.5557 
0.5536 
0.5513 
0.5487 
0.5459 
0.5429 
0.5398 
0.5366 
0.5334 
0.5302 
0.5271 
0.5241 
0.5216 
0.5199 
0.5181 
0.5163 
0.5145 
0.5127 
0.511 
0.5094 
0.5077 
0.5068 
0.5068 
0..5067 
0.5065 
0.5064 
0.5063 
0.5059 
0.5054 
0.5049 
0.5044 
0.5039 
0.5038 
0.5042 
0.5045 
0.5049 
0.5055 
0.5062 
0.507 
0.5079 
0.509 
0.5102 
0.5115 



0.004165 
0.005846 
0.008194 
0.01146 

0.016 
0.02228 
0.03093 
0.04278 
0.05893 
0.08081 
0.1102 
0.1496 
0.2017 
0.2702 
0.3596 
0.4744 
0.6178 
0.8018 
1.036 
1.331 

I. 702 
2.166 
2.741 
3.453 
4.303 
5.306 
6.522 
7.986 
9.724 

II. 8 
14.27 
17.18 

20.6 
24.57 
29.15 
34.37 
40.27 
46.99 
54.61 
63.25 
72.97 
83.87 
96.04 
109.6 
124.6 
141.3 
159.6 
179.9 
202.3 



0.1659 
0.1856 
0.2076 
0.232 
0.2589 
0.2887 
0.3215 
0.3575 
0.3968 
0.4395 
0.4858 
0.53-56 
0.5888 
0.64.56 
0.70-56 
0.7684 
0.8323 
0.8996 
0.9699 
1.043 
1.119 
1.197 
1.277 
1.36 
1.442 
1.522 
1.603 
1.685 
1.766 
1.848 
1.929 
2.01 
2.09 
2.168 
2.244 
2.316 
2.384 
2.449 
2.511 
2.569 
2.624 
2.675 
2.722 
2.765 
2.804 
2.839 
2.87 
2.897 
2.92 



0.002009 
0.002816 
0.003941 
0.005504 
0.00767 
0.01066 
0.01477 
0.02039 
0.02805 
0.03842 
0.05239 
0.07111 
0.09605 
0.1291 
0.1725 
0.2294 
0.3033 
0.399 
0.5219 
0.679 
0.8779 
1.128 
1.443 
1.837 
2.327 
2.934 
3.682 
4.6 
5.72 
7.082 
8.729 
10.71 
13.09 
15.93 
19.3 
23.28 
27.97 
33.46 
39.88 
47.37 
56.06 
66.14 
77.81 
91.29 
106.8 
124.7 
145.2 
168.7 
195.5 



0.1302 
0.1456 
0.1628 
0.1819 
0.203 
0.2263 
0.252 
0.2802 
0.311 
0.3447 
0.3814 
0.4211 
0.464 
0.5101 
0.5-595 
0.6123 
0.6684 
0.7278 
0.7906 
0.8565 
0.9251 
0.996 

1.07 
1.146 
1.225 
1.306 

1.39 
1.475 
1.562 
1.651 

1.74 
1.831 
1.923 
2.015 
2.108 

2.2 
2.293 
2.385 
2.477 
2.569 

2.66 
2.752 
2.842 
2.933 
3.024 
3.114 
3.205 
3.295 
3.386 



0.5972 
0.5967 
0.5961 
0.5954 
0.5946 
0.5936 
0.5924 

0.591 
0.5894 
0.5876 
0.5856 
0.5833 
0.5809 
0.5782 
0.5753 
0.5723 
0.5692 

0..566 
0.5627 
0.5594 
0.5561 

0.553 

0.55 

0.547 
0.5441 
0.5412 
0.5383 
0.5355 
0.5327 
0.5299 

0.527 
0.5242 
0.5213 
0.5183 
0.5153 
0.5121 
0.5089 
0.5054 
0.5018 
0.4981 
0.4942 
0.4901 
0.4859 
0.4816 
0.4772 
0.4727 

0.468 
0.4633 
0.4584 
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4.4. Tabulating the Ternary Diagram 

For the 3-layer model of solid exoplanet, points of a curve segment on the ternary 
diagram represent all the solutions for a given mass-radius input. These ternary diagrams 
are tabulated (Table [3]) with the intent to make comparison to observations easier. 

Usually, there are infinite combinations (solutions) of Fe, MgSiOs and H2O mass frac- 
tions which can give the same mass-radius pair. All the combinations togeth er form a curve 



segment on the ternary di agram of Fe, MgSiOs and H2O mass fractions (jZeng fc Seagei 



20081 : [Valencia et al.ll2007af ). This curve segment can be approximated by 3 points on it: two 
endpoints where one or more out of the 3 layers are absent and one point in between where 
all 3 layers are present to give the same mass and radius. The two endpoints correspond 
to the minimum central pressure (pOmm) and maximum central pressure [pOmax) allowed 
for the given mass-radius pair. The middle point is chosen to have the central pressure 

Table [3] contains 8 columns: 

1st column: Mass. The masses range from 0.1 to 100 with 41 points in total. 
The range between 0.1 and 1 is equally divided into 10 sections in logarithmic scale. 
The range between 1 and 10 is equally divided into 20 sections in logarithmic scale. And 
the range between 10 and 100 is equally divided into 10 sections in logarithmic scale. 

2nd column: Radius. For each mass M in Table [3] there are 12 radius values, 11 
of which are equally spaced within the allowed range: RpeiM) + {Rh2o{M) — RpeiM)) * i, 
where i = 0, 0.1, 0.2, 0.9, 1.0. The 12-th radius value (RMgSiOsiM)) is inserted into the list 
corresponding to the pure-MgSiOs planet radius (see Table [1]) for mass M. Here Rpe^M), 
RhigSiOsiM), and Rh2o{M) are the radii for planets with mass M composed of pure-Fe, 
pure-MgSiOs, and pure-H20 correspondingly. 

Overall, there are 41 * 12 = 492 different mass-radius pairs in Table [31 For each (M, R), 
3 cases: pO^m, pOmid, and pO^a^; are listed. 

3rd column: central pressure pO (Pascal) in logarithmic base- 10 scale. 

4th column: pl/pO, the ratio of pi (the first boundary pressure, i.e., the pressure at the 
Fe-MgSiOs boundary) over pO. 

5th column: p2/pl, the ratio of p2 (the second boundary pressure, i.e., the pressure at 
the MgSi03-H20 boundary) over pi. 

6th column: Fe mass fraction (the ratio of the Fe-layer mass over the total mass of the 
planet). 
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7th column: MgSiOs mass fraction (the ratio of the MgSiOs-layer mass over the total 
mass of the planet). 

8th column: H2O mass fraction (the ratio of the H20-layer mass of over the total mass 
of the planet). 

6th, 7th and 8th columns always add up to one. 



Table 3: Table for Ternary Diagram 



M(Me) 


R(i?e) 


logio(pO) 


pl/pO 


p2/pl 


Fe 


MgSiOa 


H2O 


0.1 


0.3888 


10.9212 








1 








0.1 


0.3888 


10.9212 








1 








0.1 


0.3888 


10.9212 








1 








0.1 


0.4226 


10.9066 


0.133 





0.759 


0.241 





0.1 


0.4226 


10.9126 


0.102 


0.046 


0.818 


0.172 


0.01 


0.1 


0.4226 


10.9186 


0.024 


1 


0.953 





0.047 


100 


4.102 


13.0979 


1 


1 








1 



A dynamic and interactive tool to characterize and illustrate the interior structure of 
exoplanets built upon Table [3] and other models in this paper is available on the website 
|http : / / www . cf a . harvard . edu/ ~ Izeng, 



4.5. Generate curve segment on ternary diagram using Table [3] 



One utility of Table Ois to generate the curve segment on the 3-layer ternary diagram for 
a given mass-radius pair. As an example, for M=1M® and R=1.0281i?g3, the table provides 
3 pO's. For each pO, the mass fractions of Fe, MgSiOs, and H2O are given to determine a 
point on the ternary diagram. Then, a parabolic fit (see Fig. [5]) through the 3 points is a 
good approximation to the ac tual curve segment. T his parabola may intersect the maximum 
collisional stripping curve by Marcus et al.l ( j2010l ). indicating that the portion of parabola 
beneath the intersection point may be ruled out by planet formation theory. 
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Mass:1M^ Haclius:1.0281R 




100 



Fig. 5. — The red, orange and purple points correspond to logio{pO {in Pa)) = 11.4391, 
11.5863, and 11.7336. The mass fractions are (1) red point: FeMF=0.083, MgSiO3MF=0.917, 
H2OMF=0; (2) orange point: FeMF=0.244, MgSiO3MF=0.711, H2OMF=0.046; (3) purple 
point: FeMF=0.728, MgSiO3MF=0, H2OMF=0.272. MF here stands for mass fraction. The 
blue cur ve is the parabolic fi t. The red dashed curve is the maximum collisional stripping 
curve by iMarcus et al.l ( 120101 ) . 
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5. Conclusion 

The 2-layer and 3-layer models for solid exoplanets composed of Fe, MgSiOa, and H2O 
are the focus of this paper. The mass-radius contours (Fig. [2]) are provided for the 2-layer 
model, useful for readers to quickly calculate the interior structure of a solid exoplanet. The 
2-parameter contour mesh may also help one build physical insights into the solid exoplanet 
interior structure. 

The complete 3-layer mass-fraction ternary diagram is tabulated (Table [3]), useful for 
readers to interpolate and calculate all solutions as the mass fractions of the 3 layers for a 
given mass-radius input. The details of the EOS of Fe, MgSiOa, and H2O and how they are 
calculated and used in this paper are discussed in section [3] and shown in Fig. [TJ 

A dynamic and interactive tool to characterize and illustrate the interior structure of 
exoplanets built upon Table [3] and other models in this paper is available on the website 
http : / / www . cf a . harvard . edu/ ~ Izeng 

The effect of Fe partitioning between mantle and core on mass-radius relation is explored 
in section 14. 3[ and the result is shown in Fig. H] and Table [2l 

With the ongoing Kepler Mission and many other exoplanet searching projects, we hope 
this paper could provide a handy tool for observers to fast characterize the interior structure 
of exoplanets already discovered or to be discovered, and further our understanding of those 
worlds beyond our own. 
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